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ABSTRACT 

I 

Core  physics  measurements  were  performed  in  November  1969  on  the 
first  Lefueled  core  of  the  MH-lA  (Sturgis).  Th4$  report  presents  and 
anal;,  -  .»  the  data  recorded.  The  basic  physical  parameters  describing 
the  cere  are  derived  and  compared  with  the  pre-refueling  values  and 
the  theoretical  values. 
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INTRODUCTION 


X 

Core  physics  measurements  were  performed  on  the  refueled 
MH-1A  core  in  November  1969.  The  results  are  analyzed  herein. 
Those  parameters  which  are  important  in  understanding  cere 
behavior  (such  as  the  temperature  and  power  coefficients  of 
reactivity,  rod  worths,  xenon  build-up 4  critical  bank  position) 
are  deiived  from  the  data.  Comparisons  are  made  with  expected 
valut  or  these  quantities. 


II,  DISCUSSION 

Ii  .  ial  Approach  to  Criticality 

The  initial  approach  to  criticality  after  the  first  refueling 
of  the  MH-1A  (STURGIS)  was  conducted  on  11  November  1969.  Table  1 
gives  the  data  for  the  approach  and  Figure  I  gives  the  inverse 
multiplication  curve.  The  curves  for  both  channels  predicted  a 
critical  bank  of  11.5*"  inches  withdrawal.  The  measured  critical 
bank  was  11.48  inches.  This  was  at  a  temperature  of  150°F  and  a 
pressure  of  320  psig. 

Turing  a  two-decade  power  rise  after  achieving  criticality, a 
period  of  87.3  seconds  for  a  12  rod  bank  of  11.53  inches  withdrawal 
was  observed.  This  corresponds  to  a  positive  reactivity  of  10.3  c 
or  nb~ut  $2.06  per  inch  differential  bank  worth  at  11.5  inches. 

The  rod  latch  check  was  performed  using  the  reactivity 
computer.  All  twelve  rods  were  found  to  be  latched.  Verification 
of  rod  movement  required  less  than  .5  inch  insertion  of  each  rod 
to  produce  an  observable  negative  reactivity. 

TABLE  I 


Initial  Approach  to  Criticality 


12  Rod  Bank 
Position  (inches) 

Source  Range 
Channel  1  (cpm) 

Source  Range 
Channel  2  (cpm) 

0.03 

1245 

1237 

2.00 

1297 

1231 

4.00 

1361 

1356 

6.00 

1516 

1429 

8.00 

2102 

2142 
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TABLE  I  (continued) 


Initial  Approach  to  Criticality 


12  Rod  Bank 
Position  (inches) 

Source  Range 

Channel  1  (cpm) 

Source  Range 
Channel  2  (cpm) 

9.00 

3031 

2834 

10.  UO 

5087 

4660 

10.50 

7398 

7289 

11.00 

14496 

13715 

B.  Stuck  Rod  Shutdown  Margin  Test 

The  stuck  rod  shutdown  margin  te9t  was  performed  on  11  November  1969. 
Two  inner  rods  and  three  out nr  rods  were  evaluated  in  this  test  (Table  II) 
The  critical  eleven  rod  bank  positions  for  the  simulated  stuck  rods 
Indicate  the  inner  rods  have  similar  shutdown  margins,  and  the  outer  rods 
have  similar  shutdown  margins.  It  was  therefore  not  necessary  to  measure 
the  shutdown  margin  for  each  rod  but  only  for  representative  ones. 

The  critical  eleven  rod  bank  for  an  outer  stuck  rod  was  slightly 
higher  than  that  for  an  inner  rod.  The  shutdown  margin  measurements 
for  an  inner  rod  (#1)  and  an  outer  rod  (# 5)  indicated  however,  that 
an  inner  rod  had  a  significantly  greater  margin  than  did  the  outer  rod. 

See  Figure  2  for  the  core  layout.  Pince  the  eleven  rod  critical 
positions  were  different  by  less  than  three  percent  in  the  most  extreme 
case,  the  variation  in  shutdown  margin  was  improbable  if  not  impossible. 
However,  as  a  result  of  the  extreme  space  dependence  observed  during 
the  refueling  (Ref  5)  It  was  decided  to  perform  a  second  shutdown 
margin  measurement  on  an  outer  rod  (#12)  which  was  further  from  the 
detector  supplying  the  reactivity  computer.  As  may  be  seen  from 
Table  II,  this  additional  test  further  demonstrated  the  space  dependence 
of  the  MH-1A.  Looking  at  Figure  2  one  can  see  that  as  the  "stuck  rod" 
gets  closer  to  the  detector,  the  "shutdown  margin"  decreases.  This 
is  in  direct  contradiction  to  the  critical  bank  positions. 

In  all  cases  the  required  shutdown  margin  of  -1.0  percent  Ak/k 
(-$1.37)  was  far  exceeded.  The  shutdown  margin  measured  for  rod  #5 
(-$2.60)  is  the  least  conservative  value,  and  it  is  almost  double  the 
required  margin.  It  is  also,  on  the  basis  of  engineering  judgement, 
the  most  realistic. 


TABLE  II 

STUCK  ROD  SHUTDOWN  MARGIN  MEASUREMENTS 


'iitroJ 

Rod  // 

Position 

(inches) 

11  Rod  Bank 

(inches) 

Shutdown 

Check 

Prim 

Temp  (°F) 

Prim 

Pres 

(psift) 

2 

35.84 

9.89 

152 

319 

I 

35.73 

9.89 

■m 

150 

330 

H 

35.89 

10.18 

145 

345 

5 

36.10 

10.10 

-$2.60 

147 

340 

12 

35.96 

10.00 

-$3.60 

147 

320 

C.  Control  Rod  Calibrations 


The  calibrations  of  control  rods  #1  and  #12  were  performed  on 
11  November  1969  and  15  November  1969;  the  "cold"  calibrations  being 
performed  on  the  former  date  and  the  "hot"  calibrations  on  the  latter 
date.  Each  rod  was  calibrated  singly  against  the  11- rod  bank  of  the 
remaining  rods.  The  rod  being  calibrated  was  withdrawn  an  amount 
sufficient  to  produce  a  desired  positive  reactivity  insertion.  The 
magnitude  of  the  insertion  was  then  read  from  the  reactivity 
corapt.  j  .  After  the  reading  was  completed,  the  11-rod  bank  was  moved 
so  as  to  make  the  reactor  slightly  subcritical,  and  the  negative 
reactivity  insertion  was  determined  by  the  computer.  This  process 
was  rept  :ted  until  the  rod  being  calibrated  was  fully  withdrawn. 

T«j:  Ill  compares  the  fully  withdrawn  integral  worths  for  rod  #1 

ro«J  tfl2  for  this  core  and  *ts  predecessor  (Ref  1).  Note  that  the 
information  from  the  October  1968  report  has  been  extrapolated  slightly 
so  that  it  represents  the  expected  fully  withdrawn  value.  It  should 
also  be  noted  that  since  rod  #5  and  #12  are  both  outer  rods,  it  is 
valid  for  them  to  be  compared  directly. 

The  rod  worths  measured  in  November  1969  are  vastly  different  from 
those  of  the  October  1968  report.  At  first  glance  it  would  appear  that 
it  is  impossible  to  make  any  comparisons.  This  is  not  the  case  however. 
The  reactivity  worth  ratit:  (hot  and  cold)  of  rods  #1  and  #12  are 
dii.*-tly  comparable  to  those  for  rods  #1  and  #5.  The  ratios  for  the  hot 
conditions  differ  by  less  thai  6  percent,  while  those  for  the  cold  con¬ 
ditions  differ  by  less  than  2  percent. 

'Ihis  excellent  agreement  indicates  that  the  method  which  produced 
tin  data  for  the  October  1968  report  and  the  November  1969  data  were 
consistent. 

Since  the  cores  investigated  for  each  of  the  tests  are  different, 
it  is  difficult  to  draw  any  conclusions  or  comparisons  between  actual 
numerical  values.  A  general  comment  may  be  made  however.  The  shuffled 
core  exhibits  a  much  smaller  loss  in  rod  worth  in  going  from  the  "cold" 
to  the  "hot"  condition  than  did  the  original  core  at  BOL.  The  data 
indicates  that  this  is  a  true  effect  and  is  not  due  to  a  measurement 
error 

Tables  IV  through  VII  summarize  the  measured  values  of  reactivity 
(differential  and  integral)  as  a  function  of  rod  position.  Figures  3 
and  4  indicate  the  values  of  integral  rod  worth  as  a  function  of  rod 
position. 
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TABLE  III 


FULLY  WITHDRAWN  INTEGRAL  ROD  WORTHS 


Rod  #1 

Rod  #5 

Rod  #12 

Nov  1969 

$2.84  (35.7  in,  :F4°F) 

—  —  — 

$1.53  (36.0  in,  489°F) 

(Ref  1) 

Oct  1968 

$1.92  (35.7  in.  479eF) 

$1.15  (36.0  in.  4 79 °F) 

Nov  1969 

$3.18  (35.7  in,  147°F) 

-  —  — 

$2.33  (36.0  in,  147°F) 

(Ref  1) 

Oct  1968 

$3.70  (35.7  in,  1016F) 

$2.78  (36.0  in,  101° F) 

— 

Figure  3a  and  3b  shows  the  differential  rod  worth  increases  during  the 
last  9  inches  of  travel.  This  is  explained  by  noting  that  the  last  9  inches 
of  the  control  rod  follower  have  no  burnable  poisons,  and  that  this  section 
of  the  follower  is  moving  into  a  relatively  high  flux  peak.  In  figure  4a 
this  effect  occurs  only  to  a  small  decree  since  the  bank  is  much  lower  in 
this  case.  These  two  types  of  behavior  were  also  observed  in  the  core 
physics  measurements  of  the  previous  core  (Ref  1). 
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REACTIVITY  (f/?) 


FIGURE  3.  CONTROL  ROD# I  INTEGRAL  WORTH  CURVES 
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«*Tv»  'i**1*  ‘  ■  '**»*’*«  ‘ 


(INCHES  WITHDRAWAL ) 

FIGURE  3a.  CONTROL  ROD# I  DIFFERENTIAL  WORTH  CURVE 
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T* 4«5*F 
P«I355  PSJG 


AVERAGE  ROD  POSITION 
(INCHES  WITHDRAWAL) 


FIGURE  3b.  CONTROL  ROD  #  I  DIFFERENTIAL  WORTH  CURVE 
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TABLE  IV 


Control  Rod  SI  Calibration  at  147eF,  320  psig,  11  Nov  1969 


+ 


Pod  Pos 

Ltion  (Inc 

ties  Withdrawn)  Reactivity 

| 

Before 

After 

Average 

($) 

Integral 

<s> 

c/inch 

0.05 

2.00 

1.00 

.056 

.056 

2.8 

2.00 

3.00 

2.50 

.054 

.110 

5.4 

3.00 

3.50 

3.25 

.034 

.144 

6.8 

3.50 

4.50 

4.00 

.08 

.224 

8 

4.50 

5.50 

5.00 

.10 

.324 

_  10  .. 

5.30 

6.00 

5.75 

.058 

.382 

11.6 

6.00 

7.00 

6.50 . 

.13 

.512  _ 

_  . 13  . 

7.00 

7.50 

7.25 

.07 

.582 

14 

7.50 

8.00 

7.75 

.08 

.662  . 

16 

8.00 

8.50 

8.25 

.08 

.742 

16 

8.50 

9.00 

8.75 

.08 

.822 

16 

9.00 

9.50 

9.25 

.  08 

.902 

16 

9.50 

10.00 

9.75 

.08 

.982  . 

16 

10.00 

10.50 

10.25 

.08 

1.062  . 

16 

10.50 

11.00 

10.75 

.07 

1.132 

14 

11.00 

11.50 

11.25 

.09 

18 

11.50 

12.00 

11.75 

.08 

1.302 

16 

12.00 

12.50 

12.25 

.07 

1.372 

14 

12.50 

13.00 

17.75 

.08 

1.452 

16 

13.00 

13.50 

13.25 

.08 

1.532 

16 

13.50 

14.00 

13.75 

.07 

14 

14.00 

14.50 

14.25 

.08 

1.682 

16 

14.50 

15.00 

14.75 

.06 

1.742 

12 

15.00 

15.50 

15.25 

.06 

1.802 

12 
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TABLE  IV 


Control  R,.d  /•'!  Calibration  at  147°K,  320  psig,  LI  Nov  i960 


Red  Position  (Inches  Withdrawn 


15.50 


Rtactivit 


16.00 


16.50 


17.00 


17.50 


18.00 


33.00 


34.00 


35.00 


TABLE  V 


Control  Rod  //I  Calibration  at  485°F,  1355  PSIG,  15  Nov.  1969 


Rod  Position  (Inches  Withdrawn _ _ Reactivity 


Before 

After 

_ 

Average 

($) 

$/inch 

16.25 

16.75 

16.50 

.060 

1.720 

.120 

16.75 

17.25 

17.00 

I  1 

1.775 

.110 

17.25 

17.75 

17.50 

in 

1.830 

.110 

17.75 

18.25 

18.00 

■ffll 

1.880 

.100 

18.25 

18.75 

18.50 

.05 

1.930 

.100 

18.75 

19.25 

19.00 

■n 

1.985 

.110 

19.25 

20.00 

19.63 

1 

2.050 

.087 

20.00 

21.00 

20.50 

.085 

2.135 

.085 

21.00 

22.00 

21.50 

.065 

2.200 

.065 

22.00 

23.00 

22.50 

HR 

2.255 

.055 

23.00 

24.00 

23.50 

mm 

2.295 

.040 

24.00 

25.50 

24.75 

2.345 

.033 

25.50 

27.50 

26.50 

mm 

.030 

27.50 

29.50 

28.50 

.  2.470 

.0325 

29.50 

31.00 

30.25 

2.545 

.050 

31.00 

32.00 

31.50 

.065 

2.610 

.065 

32.00 

33.00 

32.50 

.060 

2.670 

.060 

33.00 

34.00 

33.50 

.06 

2.730 

.060 

34.00 

35.00 

34.50 

.06 

2.790 

.060 

35.00 

35.68 

35.34 

.05 

2.840 

— 
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TABLE  VI 


* 


Control  Rod  #12  Calibration  at  1 4 7 ° K ,  320  PS1C,  11  Nov.  1969 


TAI 

Control  Rod  #12  Calibration  , 


j  Rod  Positiondnchea  Withdrawn) 

Before 

After 

Average 

25.00 

26.50 

25.75 

26.50 

28.00 

27.25 

28.00 

30.00 

29.00 

30.00 

■m 

31.00 

R99 

■H 

33.00 

34.00 


35.95 


34.98 


VII 


at  489°F,  1360  PSIG,  15  Nov.  1969 


Reactivity 

($) 

B9 

$/inch 

11 

1.305 

.030 

1  .050 

1.355 

• 

o 

U) 

U) 

■H 

1.410 

.028 

HH 

1.445 

.018 

mm 

1.485 

.020 

.045 

1.530 

.023 
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D.  Temperature  Coefficient  of  Reactivity 

The  temperature  coefficient  measurement  w~  performed  on  11  Nove>aber 
1969  (lower  temperature  range)  and  15  November  9  (upper  temperature 
range).  The  temperature  coefficient  was  measurec  by  establishing  a 
primary  system  linear  heat-up  rate  of  ,v60,F  per  hour  and  recording  the 
changes  in  temperature  and  reactivity  a.-.  function  of  time  during  the 
heat-up.  The  heat-up  was  accomplished  by  osing  both  pump  heat  an t  reactor 
heat.  From  this  data,  an  expression  for  both  the  inactivity  and  the 
temperature  as  a  quadratic  function  of  time  was  found,  fitting  the  data 
by  the  method  of  moments.  The  resulting  expressions  were: 

1 


R(t)  *  4. 0x10“ V  +  1.05xl0“2t  -3.5xlO~2 


T(t) 


■>-4.2 


-1. 


xlO  t  +  7.29x10  xt  +1.83x10 


11  Nov  data 


R(t)  =  2.0xl0'5t2  +  1.95xl0‘2t  -  4. 3x10“ 3 
■:(t)  =-3.7xl0“V  +  1.066t  +  2.42xl02 


15  Nov  data 


where  time  t  is  in  minutes,  reactivity  R  is  in  $,  and  temperature  T  is 
in  °F. 


From  these  fitted  curves,  the  temperature  coefficient  is  calculated 


thus : 


A  R 
/  R  t 


(T)  = 


KUft1)  -  R(t) 


T  >'  c  +  t ' )  -  T(t) 


where  t'  is  the  increment  between  successive  time  steps  (t'  =  10  minutes 
was  used)  and  T  is  the  average  temperature  over  the  time  step,  T  ■*  T(t4.5t') 

Figure  5  illustrates  the  results  Note  that  the  data  on  the  two 
days  demonstrates  contradictory  trends:  The  data  of  11  November  results 
in  a  curve  whose  slope  decreases  as  temperature  increases,  while  the 
data  of  15  November  shows  an  opposite  trend.  The  curve  of  11  November 
is  in  error  due  to  variations  in  heat-up  rate,  filling  of  the  steam 
generator,  and  lack  of  overall  experience  in  running  the  temperature 
coefficient  measurements  with  the  reactivity  computer.  The  data  of 
15  November  shows  much  better  results.  All  the  data  was  used  however, 
in  order  to  provide  information  over  the  entire  operating  range  of  the 
reactor. 


Using  the  data  of  Figure  5,  a  single  quadratic  curve  was  found, 
again  by  the  method  of  moments,  to  give  the  best  fit  to  the  data: 


Ae_ 


A  T 


=  4. 39xlO-8T2  +  2. 84x10“ ;>T  +  9.0xl0"3  ($/°F) 
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This  curve  is  plotted  in  Figure  6,  At  68° F,  it  yields  a  tenperature 
coefficient  of  reactivity  of  -l.llc/°F  or  -8.10x10“5/oj>.  At  490#F,  it 
yields  a  temperature  coefficient  of  -3.35c/°F  -2.45xlO"4/“F.  The 

temperature  defect  is  found  by  Integration  of  the  quadratic  expression: 


490 

A  p  (AT)  ■/  -£-P—  dt  *  $8.86  or  6.472  Ap 

68  AT 

Table  VIII  compares  these  results  with  previous  values. 

The  accuracy  of  fit  cf  the  quadratic  expression  to  the  data  may 
be  measured  b,  adding  up  the  reactivity  change  contributions  (see 
Tables  IX  and  X)  over  the  entire  temperature  range ,  and  comparing  this 
to  the  $8.86  obtained  above.  Since  the  data  only  begins  at  185°F, 
the  contribution  from  68°F  to  185° F  was  assumed  to  be  well  represented 
by  the  integration  of  the  quadratic  expression: 


185 

/  (.LiL)  dT  -  $1.56 

68  lT 

From  185°F  to  490° F,  the  !  2  values  were  added  up  (see  Tables  IX  and  X), 
yielding  $7.41.  Thence  the  temperature  defect  was  found  by  adding 
up  the  raw  data  result  and  the  correction  back  to  68°F. 


Ap  (AT)  -  $7.41  +  $1.56  •  $8.97 


This  result  agrees  to  within  1.2  percent  of  the  result  obtained  by 
integrating  the  quadratic  fit  (which  smooths  out  fluctuations  in  the 
data):  thus,  the  quadratic  curve  represents  well  the  true  temperature 
coefficient  behavior  in  a  smoothed-out  sense. 

TABLE  VIII 


TEMPERATURE  COEFFICIENT 

Temperature 

68°F 
490°F 

TEMPERATURE  DEFECT 

6.47%  Ap  6.4%  Ap  4.8%  Ap 


Experimental  Theoretical  (Ref  2) 

Nov  69  Oct  68  (Ref  1) 


-8.10xl0"5 

-2.45xl0“4 


-3. 0x10“ 5 
-2.95xl0"4 


-3.3x10"5/°F 
-2. 7x10”4/°F 
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COEFFICIENT  DATA 


The  data  yields  a  ip/fiT  much  larger  at  low  temperature  than  the 
theoretical  value.  However,  better  agreement  (about  16  percent  low)  is 
obtained  at  the  upper  temperature.  Though  the  value  of  the  temperature 
defect  differs  significantly  from  the  theoretical  value  (4.8  percent 
A  p) ,  it  agrees  well  with  the  value  measured  for  the  previous  core 
(6.8  percent  Ac).  Hence,  the  Ap/AT  versus  T  behavior  of  the  two 
cores  is  nearly  the  same  in  an  integral  sense,  even  though  there 
appears  to  be  a  variance  in  low  temperature  behavior. 


TABLE  IX 

temperature  Coefficient 

Date:  11  Nov  1969 


Stopwatch 

Time 

(Min) 

Average 

Primary 

Temperature 

(°F) 

Primary 
Pressure 
(psig)  _ 

Average 

12  Rod  Bank 
Position 
(Inches) 

Reactivity 

Change 

($) 

0 

185 

330 

11.92 

0 

3 

188 

330 

11.95 

.04 

8 

191 

330 

11.98 

.04 

14 

194 

330 

12.03 

.07 

19 

197 

330 

12.07 

.055 

26 

202 

33  j 

12.12 

.055 

30 

204 

335 

12.16 

.035 

37 

211 

335 

12.24 

.105 

45 

217 

335 

12.30 

,06 

50 

220 

335 

12.36 

.08 

56 

225 

335 

12.42 

.08 

61 

230 

335 

12.49 

.08 

65.5 

234 

338 

12.55 

.11 

69 

237 

340 

12.60 

.065 

75 

245 

335 

12.70 

.125 

78-80 

246 

340 

12.70 

-  .01 

83 

248 

340 

12.78 

.08 

87 

253 

340 

12.84 

.095 

** 


**  Steam  Generator  Filling 
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TABLE  IX  (CONT'D) 
Temperature  Coefficient 


Date:  11  \7>v  1‘».9 


Stopwatch 

Tine 

(Min) 

Average 

Primary 

Temperature 

(°F) 

Primary 

Pressure 

(nsif>) 

Average 

12  Rod  Bank 
Position 
(Inches) 

React ivi tv 
Change 
($) 

....  93 

259 

340 

12.93 

.10 

100 

264 

340 

13.02 

.  143 

107 

272 

350 

13.18 

•  1 

114 

279 

355 

13.33 

.183 

_120  . 

287 

370 

13.42 

.123 

128 

296 

410 

13.57 

.175 

mm 

306 

470 

13.76 

.235 

320 

540 

14.01 

.285 

160 

325 

580 

14.11 

.12 

166 

331 

610 

14.22 

.125 

172 

336 

655 

14.33 

.12 

178 

342 

690 

14.47 

.125 

184 

347 

■■■ 

14.60 

.14 

189 

350 

760 

14.71 

.125 

TABLE  X 

Temperature  Coefficient 


Stopwatch 

Time 

(Min) 


Average 

Primary 

Temperature 

(°F) 

Primary 

Pressure 

(psig) 

247 

340 

252 

340 

257 

350 

264 

380 

267 

400 

275 

400 

282 

400 

287 

400 

290 

410 

294 

435 

301 

460 

307 

480 

315 

519 

320 

549 

328 

575 

336 

575 

341 

590 

347 

630 

351 

660 

358 

700 

363 

730 

366 

760 

371 

800 

Average 
12  Rod  Bank 
Position 
inches 


12.47 


12.54 


12.62 


12.74 


12.84 


12.92 


13.04 


3, 


13.21 


13.30 


13.40 


13.51 


13.64 


13.77 


13.94 


14.08 


14.20 


14.32 


14.41 


14.56 


14.66 


14.78 


14.88 


Date:  15  Kov  1969 


Reactivity 

Change 

($> 
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TABLE  X 

Temperature  Coefficient 


Stopwatch 

Time 

(Min) 

Average 

Primary 

Temperature 

(*F) 

Primary 

Pressure 

(psig) 

Average 

12  Rod  Bank 
Position 
(inches) 

Reactivity  I 
Change 
($) 

130 

376 

830 

15.00 

.125 

135 

381 

870 

15.13 

.115 

140 

385 

900 

IS.  24 

.115 

146 

390 

950 

15.40 

.155 

152 

397 

990 

15.53 

.13 

157 

401 

1030 

15.65 

.11 

162 

405 

1080 

15.77 

.115 

166 

409 

1110 

15.86 

.085 

171 

413 

1160 

15.99 

.105 

175 

417 

1200 

16.11 

.16 

180 

421 

1240 

16.19 

.03 

185 

425 

1290 

16.33  . 

.12 

190 

429 

1335 

16.47 

.13 

196 

435 

1380 

16.63 

.15 

200 

439 

1375 

16.76 

.115 

203 

442 

1385 

16.87 

.085 

208 

447 

1385 

17.03 

.14 

212 

451 

1382 

17.20 

.13 

216 

454 

1365 

17.35 

.125 

220 

458 

1355 

17.52 

.125 

224 

462 

1365 

17.63 

.10 

227 

464 

1363 

17.78 

.13 

230 

468 

1360 

17.94 

.12 

237 

475 

1360 

18.23 

.20 
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TABLE  X 

Temperature  Coefficient 


Stopwatch 

Time 

(Min) 

Average 

Primary 

Temperature 

(°F) 

Primary 

Pressure 

(psig) 

Average 

12  Rod 
Position 
(inches) 

Reactivity 

Change 

($) 

241 

478 

1370 

18.40 

.13 

245 

482 

1358 

18.61 

.15  i 

-  -  .  -  --  j 

248 

485 

1350 

18.78 

1 

.14  | 

254 

491 

1372 

19.03 

; 

.165  i 

E.  Power  Coefficient  of  Reactivity 

The  power  coefficient  measurements  were  performed  on  19  November  1969 
(Section  I  data,  increments  increasing  power  at  minimum  xenon  conditions) 
and  on  25  November  1969  (Section  II  data,  increments  decreasing  power  from 
equilibrium  xenon  conditions) .  The  technique  employed  was  to  increase 
(decrease)  the  reactor  power  by  increasing  (decreasing)  the  generator 
load  and  then  recording  the  resulting  temperature  change.  Knowing  the 
temperature  coefficient  of  reacti.it.-,  it  is  possible  to  ascertain  the 
reactivity  change  caused  by  the  'hange  in  power.  Thus,  the  power 
coefficient  is  simply  the  ratio  of  the  reactivity  change  to  the  power 
change,  assigned  to  the  average  power  in  the  interval  of  power  change. 

Table  XI  shows  representative  recorded  information  for  these  measurements. 

The  resulting  values  of  the  power  coefficient  are  presented  in  Figure  7. 
Considerable  scatter  is  evident  in  the  calculated  values,  particularly 
those  due  to  the  Section  I  data.  This  is  due  to  the  build-up  of  xenon 
during  the  measurements,  resulting  when  a  large  amount  of  time  (4-1/2 
hours)  was  spent  in  making  the  power  changes.  Thus,  the  requirement 
of  minimum  xenon  vas  not  met  for  the  later  data  of  Section  I,  and 
these  points  (power  >  30  MW)  are  obviously  in  error  and  were  not  used  in 
computing  the  power  coefficient.  The  remaining  data  from  both  Sections 
I  and  II  was  used  in  obtaining  a  linear  curve  fit,  by  the  method  of 
moments,  shown  in  Figure  7 .  Thence  the  average  value  of  the  power 
coefficient  of  reactivity  was  found  to  be  -4.03<?/MWT  or  -2.95xl0"^/MWT. 

This  is  to  be  compared  to  the  theoretical  value  of  -2. 1x10  yMWT. (Ref  2). 

In  spite  of  the  data  scatter,  all  the  experimental  values  are 
significantly  higher  than  the  theoretical  ones. 
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POWER  COEFFICIENT 
(4 /MW) 


3*  SECTION  I  DMA 
®«9CCTI0N1I  QflTA 


(MW) 

NB.  IN  LINEAR  FIT,  THE  THREE 
HIGHEST  POWER  DATA  FONTS 
OF  SECTION  I  DATA  WERE 
DISREGARDED. 

FIGURE  7.  POWER  COEFFICIENT  DATA 
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The  measurements  of  power  coefficient  gave  poor  results.  The 
analysis  of  the  power  coefficient  data  was  rendered  difficult  and  some 
somewhat  arbitrary  by  three  lapses  in  experimental  technique: 

1.  The  measurements  were  not  performed  quickly  enough,  thus 
permitting  xenon  to  build-up  and  obscure  the  data. 

2.  On  the  data  recording  sheets,  the  readings  from  the  Keithley 
pico-amneter  were  either  omitted  or  illegible,  making  it  difficult  to 
verify  the  power  changes  or  to  normalize  the  readings  of  the  power 
channu '  >> . 


J.  A  recalibration  was  performed  in  the  middle  of  the  measurements. 


TABLE  XI 

POWER  COEFFICIENT  DATA 


Time 

Primary 

Average 

Average 

Remarks 

(Hours ) 

Pressure 

Rod  Bank 

Primary 

(psig) 

Position 

(Inches) 

Temp  (°F) 

BH 

Start  of 

1300 

1375 

20.60 

518 

13 

Section  1  Data 

(19  Nov  69) 

13JO 

1367 

20.60 

515 

19.5 

1400 

1368 

20.60 

498 

54 

1 4  30 

1345 

20.60 

488 

71 

1555 

- 

- 

- 

71 

Recalibration  of 
Power  Level 

1555 

1350 

21.55 

514 

61 

1610 

1345 

21.55 

502 

75 

1633 

1365 

21.55 

492 

85 

1656 

1315 

21.55 

493 

83.5 

1728 

1370 

21.55 

483 

98 

1334 

1365 

26.70 

490 

95 

Start  of  Section  II 
Data  (25  Nov  69) 

1338 

1360 

26.70 

498 

80 

1350 

1362 

26.07 

490 

80 

1358 

1363 

26.07 

499 

60 

*  Channel  7  Data  used,  considered  the  most  consistent. 
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TABLE  *T (CONT’D) 
Power  Coefficient  Data 


Time 

(Hours) 

Primary 

Pressure 

(psig) 

Averate 

Red  Bank 
Position 
(Inches) 

Average 
Primary 
Temp  (°F) 

Power 

Level 

*00 

Remarks 

1413 

1377 

25.45 

490 

61 

1421 

1360 

25.45 

502 

42 

1435 

1374 

24.75 

490 

42 

Jl 

504 

15 

*  Channel  7  data  used,  considered  the  most  consistent. 

F.  Equilibrium  Xenon 

The  xenon  equilibrium  measurements,  with  the  reactor  at  near-full 
power,  were  performed  on  19  November  1969  through  20  November  1969. 

The  technique  employed  for  these  measurements  used  the  previously 
measured  temperature  coefficient  to  ascertain  the  amount  of  poisoning 
caused  by  xenon  build-up.  As  xenon  built-up,  the  loss  of  reactivity 
caused  by  the  increased  absorption  of  neutrons  was  balanced  by  the 
decrease  in  temperature  and  power.  Periodically  the  control  rods  were 
withd-  ti  to  return  the  mean  primary  temperature  to  ^490oF  and 
subsequently  raise  the  power  level.  By  summing  the  changes  in  primary 
temperature,  multiplying  this  sum  by  the  temperature  coefficient,  and 
adding  contributions  due  to  power  changes,  the  negative  reactivity 
attributable  to  xenon  build-up  was  calculated. 

Due  to  an  omission  in  recording  data,  the  rod  bank  was  moved 
(at  time  =  18.13  hours)  without  recording  the  temperature;  thus,  this 
datum  point  could  not  be  calculated.  Instead,  a  method  of  moments 
linear  fit  (on  semi-logrithmic  paper,  since  the  expected  exponential 
form  would  yield  a  straight  line)  was  used  to  predict  this  datum 
(see  Figure  8) .  Succeeding  points  depended  on  this  datum  since  the 
xenon  build-up  was  calculated  in  an  integral  manner. 

In  addition  to  the  error  associated  with  the  omitted  data,  is 
the  error  in  overall  loss  of  reactivity  due  to  xenon  build-up  caused 
by  the  amount  of  time  consumed  in  performing  the  power  coefficient  test. 
Since  the  power  coefficient  data  is  so  scattered,  the  total  amount  of 
reactivity  loss  due  to  xenon  build-up  during  the  power  coefficient 
measurement  is  hard  to  estimate.  In  addition,  the  time  of  equivalent 
full-power  operation  is  also  hard  to  estimate.  From  the  data  of  Table  XI, 
the  estimated  time  of  full-power  operation  can  be  calculated  by  summing 
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MD*  MISSING  DATUM 
A  VALUES  PREDICTED  FROM  GRAPH 
GIVEN  VALUE  FORT  -  ».I3 
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FIGURE  8.  EXTRAPOLATION  TO  OBTAIN  MISSING  DATUM 


* 


the  fractional  power  operation  during  each  time  interval.  This  yields 
2.95  hours  of  equivalent  full-power  operation  during  the  power 
coefficient  measurements,  equivalent  to  building  in  approximately 
22c  of  xenon.  Figure  9  shows  the  resulting  time  shifted  data.  For 
purposes  of  comparisons  the  theoretical  curve  is  also  presented  (see 
below  for  how  it  was  obtained).  It  is  apparent  that  the  data  values 
are  significantly  smaller  than  the  theoretical  values.  The  equilibrium 
value  obtained  from  the  data  is  -$1.85  or  -1.4  percent  ip  (obtained  by 
adding  22c  to  the  asymptotic  value  indicated  by  the  data  in  Figure  9). 

It  is  meaningful  to  calculate  what  the  number  density  of  Xe-135 
as  a  function  of  time  would  be  from  theoretical  considerations,  and 
compare  it  to  the  behavior  observed  in  the  data.  Now,  Xe-135  is 
formed  by  the  following  decay  chain: 


fission  fission 

Te-135  -  1-135  ■*  Xe-135  Cs-135 


Because  Te-135  decays  very  rapidly,  the  usual  approximation  (3)  is  made 
that  1-135  is  directly  produced  in  fission. 

The  theoretical  xenon  build-up  curve  may  be  obtained  by  solving 
the  following  differential  equations:  (Ref  3), 

-  AI  +  B 
dt 

4~  **  CX  +  X  I  +  D 
dt  I 

where 

A  - 

3  -  yf  f* 

C  *  —  ( A  +  a  $) 
x  ax 

D  =  y  E  <p 
x  f 

I  *  1-135  number  density 
X  =  Xe-135  number  density 
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FIGURE  9.  XENON  BUILD •  UP  (EXPERIMENTAL  DATA) 


The  resulting  expression  for  the  time-dependent  xenon  number  density  is: 


X(t) 


-(.'  +o  -*)t 
tl-e  x  ax  j 


-(>.  +0  <*-)t  ->Tt.  . 
=>  v  v  T  I  * 


The  foil. wing  parameters  were  used  (Reference  3): 

=  2. 7x10“ cm" 

,  -  0.061 
I 

V  =  0.003 
x 

A  -  0.1035  hr'1 
*  0.0753  hr-1 

And  the  fission  cross-section  and  number  density  for  uranium  and  the 
thermal  llux  were  obtained  from  the  LEOPARD  program  (Ref  4) : 

7  -N  =  (3.85xl0"4/A3)  (282xlO~24  cm2)  (108  A/cm)3  -  0.108  cm-1 

f  235  *235 

$  =  1.  j>:1033/cm  sec 

ITie  of235  u.:ed  was  the  Wigner-Wilkins  spectrum  averaged  value  given  by 
LEOPARD.  Using  these  parameters,  the  solution  for  X(t)  becomes: 

X(t)  =  0.0705xl0-15  [l-e~  32t]  +  3.35xl015  [.571  e_*182t  -e_'104t3  +  1.43xl015 

where  the  units  of  time  are  hours,  and  X  is  in  ///cm3. 

Now  the  reactivity  p  is  defined  (Reference  3) 

p  =  1ax/_£_f  =  ~Xnax/  V  f 

vpe  vpe 

From  LEOPARD,  p  =  0.75.  Also  v  =  2.44  and  e  =  1.0  were  assumed. 
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The  resulting  time  dependent  behavior  of  d  as  xenon  build«-up  is 
shown  in  Figure  10  and  labeled  "theoretical  curve." 

The  data  as  obtained  from  the  measurements  of  19  November  1969  are 
also  shewn  without  the  time  shift  on  Figure  10,  labeled  "uncorrected  data." 
It  is  clear  that  this  curve  falls  far  below  the  theoretical  one.  The 
results  are  very  sensitive  to  the  measurements  of  temperature  made,  sirce 
the  reactivity  is  calculated  in  a  summation  process. 


Thus,  it  was  assumed  that  the  temperature  drops  as  recorded  did 
not  reflect  the  true  build-up  of  xenon.  Instead,  each  degree  of  AT  was 
assumed  to  have  been  measured  to  be  an  additive  constant  5/3  less  than 
its  true  value  (e.g.  a  U  of  3  degrees  wa3  in  reality  a  £  of  (3+0  degrees). 
The  A  power  measurements  were  assumed  correct.  Using  €  as  an  unknown,  the 
total  build-up  of  xenon  was  calculated  and  set  to  equal  the  theoretical 
value,  thus  yielding  one  equation  in  the  unknown  6: 


t-40 

Z 

t*o 


t=40 

(AT*  44“  )  +  >■ 

t=o 


theoretical 


This  yields  a  value:  6  *  2.27°F. 


Thus,  the  AT  recorded  on  19  November  1969  were  modified: 

A  T'  =  AT  +  (-^4)  A 

and  the  data  corrected  accordingly.  The  resulting  curve  is  shown  in 
Figure  10,  labeled  "6-corrected  data".  It  appears  to  agree  well 
with  the  theoretical  curve. 


The  equilibrium  xenon  is  calculated  from  the  theoretical  curve  to 
be  -2.06  percent  or  -$2.82.  The  6-corrected  curve  predicts  an  equilibrium 
value  of  -$2.8. 

Ill  III.  CONCLUSIONS 


The  core  physics  tests  performed  in  November  3969  on  Core  2  of  the 
MH-.IA  have  provided  a  great  deal  of  useful  information.  They  demonstrated 
the  enormous  time  savings  which  result  from  the  use  of  a  reactivity  com¬ 
puter  in  rod  latch  verification,  stuck  rod  margin  measurement,  control 
rod  calibration,  and  temperature  coefficient  measurement.  They  also 
showed  that  the  reactor  operators  and  data  takers  must  be  meticulous 
and  concientious  in  the  performance  of  the  tests  and  in  recording  data. 
Some  of  the  data  gathered  in  these  tests  did  not  provide  the  information 
which  it  should  have.  This  is  due  to  many  factors;  poorly  worded 
procedures,  failure  to  record  required  data,  and  general  inexperience 
in  the  performance  of  core  physics  tests.  On  the  basis  of  these  tests 
the  procedural  problems  can  be  cleared  up  prior  to  the  next  tests.  The 
failure  to  record  required  data  can  be  taken  care  of  by  better  instruction 
and  orientation  of  operating  personnel.  The  problem  of  inexperience  is 
not  an  easy  one  to  overcome,  due  to  the  continual  changeover  of  plant 
and  program  personnel.  It  definitely  illustrates  the  requirement  for 
the  presence  of  an  experienced  nuclear  plant  test  engineer  in  addition 
to  the  regular  plant  operating  and  management  personnel. 
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TIME  (HR) 

FIGURE  10.  XENON  BUILD-UP  (CORRECTED  DATA) 


Table  XII  compares  the  present  results  with  those  of  previous 
measurements  and  those  expected  from  theoretical  calculations.  The 
agreement  of  present  results  with  those  based  on  previous  experience 
or  theory  appear  satisfactory. 

Even  though  some  of  the  tests  provided  poor  results,  the  entire 
test  has  resulted  in  a  better  understanding  of  the  MH-1A  and  its  Nuclear 
performance. 


TABLE  XII 

SUMMARY  COMPARISON 


Experimental 

Theoretical  (2) 

Previous  Core  (1) 

Temperature  coefficient 

(68°F) 

-0.81xl0~4/°F 

-0.33x10-4/°F 

-0. 3x10~4/°F 

(490°F) 

-2.45xlO“4/°F 

-2.7x10'4/°F 

-2.95x10~4/°F 

Power  coefficient 

-2.95x10_4/MWT 

-2. 1x10-4/MWT 

-2.01x10“4/MWT 

Temperature  defect 

6.47%  Ap 

4.8%  Ap 

6.4%  Ap 

Equilibrium  xenon 

-1.4%  Ap 

-2.1%  Ap 

-2.0%  Ap 

Critical  bank  height  (68°F) 

11.48  inches 
(@  150°F) 

11.18  inches 
(@68°F) 

11.3  inches 
(@100°F) 
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